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ABSTRACT 

We explore the constraints on those extensions to the standard models of cosmology and particle physics 
which modify the early- Universe, radiation-dominated, expansion rate S = H' /H (parametrized by the effective 
number of neutrinos N„). The constraints on S (N v ) and the baryon density parameter rjs = («b /w 7 ) = 10~ 10 ?7io, 
qq ' derived from Big Bang Nucleosynthesis (BBN, t ~ 20 minutes) are compared with those inferred from the 

Cosmic Microwave Background Anisotropy spectrum (CMB, t ~ 400 kyr) and Large Scale Structure (LSS, 
t ~ 14 Gyr). At present, BBN provides the strongest constraint on N^ (N v = 2.4 ± 0.4 at 68% confidence), but a 
CNJ ' weaker constraint on the baryon density. In contrast, while the CMB/LSS data best constrain the baryon density 

£H \ (r/io = 6.1^o j at 68% confidence), independent of N„, at present they provide a relatively weak constraint 

on N v which is, however, consistent with the standard value of N y = 3. When the best fit values and the 
allowed ranges of these CMB/LSS -derived parameters are used to calculate the BBN-predicted primordial 
abundances, there is excellent agreement with the observationally inferred abundance of deuterium and good 
agreement with 4 He, confirming the consistency between the BBN and CMB/LSS results. However, the BBN- 
predicted abundance of 7 Li is high, by a factor of 3 or more, if its observed value is uncorrected for possible 
dilution, depletion, or gravitational settling. We comment on the relation between the value of N„ and a possible 
i-C i anomaly in the matter power spectrum inferred from observations of the Ly-a forest. Comparing our BBN and 

CMB/LSS results permits us to constrain any post-BBN entropy production as well as the production of any 
q non-thermalized relativistic particles. The good agreement between our BBN and CMB/LSS results for N„ and 

5_h ■ rjs permits us to combine our constraints finding, at 95% confidence, 1.8 < N, y < 3.2 and 5.9 < 7710 < 6.4. 

Subject headings: neutrinos — early universe — expansion rate — baryon density 

d 

1 1. INTRODUCTION 

. The standard models of particle physics and of cosmology with dark energy, baryonic matter, radiation (including three species 
of light neutrinos), and dark matter is consistent with cosmological data from several widely-separated epochs. However, there 
T _j- ' is room to accommodate some models of non-standard physics within the context of this well tested, ACDM cosmology. One 
possibility, explored here, is that of a non-standard expansion rate (S = H' /H, where H is the Hubble parameter) during the 
• ' early, radiation-dominated evolution of the Universe driven perhaps, but not necessarily, by a non-standard content of relativistic 
£!D \ particles (p' R / pr), parametrized by the equivalent number of additional neutrinos (p' R = Pr + AN^p^, where AN,, = N^ - 3 prior 
to the epoch of e^ annihilation). In the epoch just prior to e^ annihilation, which is best probed by Big Bang Nucleosynthesis 
(BBN), p R = p 1 +p e +3p u = 43p 7 /8, so that 
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^ , In the epoch after the completion of e^ annihilation, best probed by the Cosmic Microwave Background (CMB) and by Large 
Scale Structure (LSS), the relations between /?r and p 1 and between S and AN,, differ from those in eq. 1, as described below in 
some detail in §2. 

We use the comparison between the predicted and observed abundances of the light elements produced during BBN, and 
between the predicted and observed CMB anisotropy spectrum, along with data from LSS observed in the present/recent Universe, 
to constrain new physics which leads to a non-standard, early-Univer se expansion rate (S ) or, equivalen tly, to place bounds on the 
effective number of neutrinos (N^); for related earlier work see, e.g., Bar ger et"ai1 (120031) : ICyburj (l2004h . In addition, the baryon 
density and, any variation in it over widely-separated epochs in the evolution of the Universe, are constrained simultaneously 
with N„, thereby testing the standard-model expectation that the ratio (by number) of baryons to CMB photons = «b/«7) 
should be unchanged from e ± annihilation (T < 1/2 MeV; t > 3 s) until the present (T = 2.725 K w 2 x 10" 10 MeV, t w 14 Gyr). 

In §2 a non-standard expansion rate (S) is related to a non-standard radiation density, parametrized by an effective number of 
neutrinos N^. In §3, the observationally inferred primordial light element (D, 4 He, 7 Li) abundances are used to constrain the 
radiation density (expansion rate) and the baryon density. In §4, we assume a flat, ACDM cosmology and use data from the 
CMB, along with a prior on the Hubble parameter from the HST key project, luminosity distances of type la supernovae, and the 
matter power spectrum to provide independent constraints on the radiation density and the baryon density. The results from these 
two epochs, very widely-separated in time, are compared constraining any post-BBN entropy production. The good agreement 
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between them permits us to combine them to obtain a joint constraint on the baryon density parameter (?7b) and the effective 
number of neutrinos (N, y ). 

2. NON-STANDARD EARLY-UNIVERSE EXPANSION RATE OR RADIATION DENSITY 

In the radiation-dominated early universe (ptot — v Pr) the expansion rate (H) is related to the radiation density through the 
Friedman equation. 



H 2 = ~—Gp R . (2) 

Any modification to the radiation density, or to the Friedman equation by a term which evolves like the radiation density (as 
the inverse fourth power of the scale factor) can be parametrized by an equivalent number of additional neutrinos AN„ where, 
prior to e ± annihilation, AN„ = N„ — 3. For the standard models of particle physics and cosmology, in the epoch after muon 
annihilation (T < 100 MeV) and prior to annihilation (T > 0.5 MeV), the radiation consists of an equilibrium mixture of 
photons, relativistic pairs, and three flavors of extremely relativistic, left-handed neutrinos (and their right-handed, antineutrino 
counterparts). In this case, the total radiation density may be written in terms of the photon density as 

43 

Pr= — p 7 = 5.375p 7 . (3) 
In this same epoch, prior to <? ± annihilation, a modified radiation density can be written as 

/°r = /°r(i + -^^ =p R (l + 0.163AN„) (4) 

where /5r is the standard-model radiation energy density and p' R is the modified, non-standard model radiation energy density. In 
a sense, this modified energy density is simply a proxy for a non-standard expansion rate during the radiation-dominated epoch 
relevant for comparison with BBN. 
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After e annihilation the surviving relativistic particles are the photons (which will redshift to the presently observed CMB) 
and the now decoupled, relic neutrinos. In the approximation that the neutrinos are fully decoupled at e ± annihilation, the 
post-e* annihilation photons are hotter than the neutrinos by a factor of T^/T v = (1 1 /4) 1 / 3 and 



PR = 

so that 
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p 7 = (l + 3x0.227)p 7 =1.681/9 7 , (6) 



p' R = p 7 (l. 681 +0.227 AN„) = p R (l +0.135AN„), (7) 

where, as before, AN„ = N^ -3. 

Ho wever, it is well known dDicus et al.ll982l:lDodelson & Turner! 1992tlHannestad & Madsenll995HLopez et al.ll999HDolg"ov1 
2002) that the standard-model neutrinos are not fully decoupled at e^ 1 annihilation. As a result, the relic neutrinos share some 
of the energy/entropy released by <? ± annihilation and they are warmer than in the fully decoupled approximation, increasing 
the ratio of the post-e ± annihilation radiation density to the photon energy density. While the post-e ^ annihilation phase sp ace 
distribution of the decoupled neutrinos is no longer that of a relativistic, Fermi-Dirac gas, according to Ma ngano et alJ (l2005h the 
additional contribution to the total energy density can be accounted for by replacing N„ = 3 with N„ = 3.046, so that 

Pr -> (1 + 3.046 x 0.227)p 7 = 1.692p 7 . (8) 

For deviations from the standard model that can be treated as equivalent to contributions from fully decoupled neutrinos, 

p^ = p R (l+0.134AN„), (9) 

where AN„ = NJ,- 3.046 in the post-e ± annihilation Universe relevant for comparison with the CMB and LSS. Note that in 
the standard model, where N„ = 3 prior to e ± annihilation, the neutrino contribution to the post-e* annihilation radiation energy 
density is equivalent to = 3.046, so that for the standard models of cosmology (standard expansion rate) and of particle physics 
(standard radiation energy density), N„ = 3, N' v = 3.046 and, AN„ = 0. 

We emphasize that although the non-standard radiation density (expansion rate) has been parametrized as if it were due to 
additional species of neutrinos, this parametrization accounts for any term in the Friedman equation whose energy density varies 
as a" 4 , where a is the scale factor. From this perspective, AN„ could either be positive or negative; the latter does not necessarily 
imply fewer than the standard-model number of neutrinos but could, for example, be a sign that the three standard-model neutrinos 
failed to be fully populated in the early Universe or, could reflect modifications to t he 3+1 dimensional Friedman equations 
arisin g from high er-dimensional extensions of the standard model of particle physics (Randall & Sundrumlll999l : lBinetruy et alJ 
2000: ICline et"al]l2000h . 
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3. AND r/ B FROM BBN 

The stage is being set for BBN when the Universe is about a tenth of a second old and the temperature is a few MeV. At this 
time the energy density of the universe is dominated by relativistic particles. When the temperature drops below ~ 2 MeV, the 
neutrinos begin to decouple from the photon-e^ 1 plasma. However, they do continue to interact with the neutrons and protons 
through the charged-current weak interactions (n+v e «-> p+e~,p+D e <-> n+e + , n <-» p+e~ + v e ), maintaining the neutron-to-proton 
ratio at its equilibrium value of n/p = exp(-Am/kT), where Am is the neutron-proton mass difference. When the temperature 
drops below ~ 0.8 MeV, and the Universe is ~ 1 second old, the reactions regulating the neutron-proton ratio become slower 
than the universal expansion rate (r„A < H). As a result, the neutron-proton ratio deviates from (exceeds) its equilibrium value, 
so that n/p > exp(-Am/kT), and the actual n/p ratio depends on the competition between the expansion rate (H) and the 
charged-current weak interaction rate (T W k), as well as on the neutron decay rate, 1 /r„, where r„ is the neutron lifetime. 

Although nuclear reactions such as n + p <-> D + 7, proceed rapidly during these epochs, the large 7-ray background (the blue- 
shifted CMB) ensures that the deuterium (D) abundance is very small, inhibiting the formation of more complex nuclei. The more 
complex nuclei begin to form only when T < 0.08 MeV, after e ± annihilation, when the Universe is about 3 minutes old. At 
this time the number density of photons with sufficient energy to photodissociate deuterium is comparable to the baryon number 
density and various two-body nuclear reactions can begin to build more complex nuclei. Note that the neutron-to-proton ratio has 
decreased slightly since "freeze-out" (at T < 0.8 MeV) through the residual two body reactions as well as via beta decay. Once 
BBN begins neutrons and proton combine to form D, 3 He, and 4 He. The absence of a stable mass-5 nuclide presents a road- 
block to the synthesis of heavier elements in the expanding, cooling Universe, ensuring that the abundances of heavier nuclides 
are severely depressed below those of the lighter nuclei. In sta ndard BBN (SB BN) only D, 3 He, 4 He, and 7 Li are produced in 
astrophysically interesting abundances (for a recent review see Steigman d2007l) ). While the BBN-predicted abundances of D, 
3 He, and 7 Li are most sensitive to the baryon density, that of 4 He is very sensitive to the neutron abundance when BBN begins 
and, therefore, to the competition between the weak-interaction rate and the universal expansion rate. The primordial abundances 
of D, 3 He, or 7 Li are baryometers, constraining r/s, while 4 He mass fraction (Yp) is a chronometer, depending mainly on 5 or, N, y . 
In the standard model of particle physics and cosmology with three species of neutrinos and their respective antineutrinos, the 
primordial element abundances depend on only one free parameter, the baryon density parameter, the post-e 111 ratio (by number) 
of baryons to photons, rjs = n%/n 1 . Thi s parameter ma y be related to J7b, the present-Universe ratio of the baryon mass density 
to the critical mass-energy density (see Steigman (2006)) 

7710= 10 10 n B /« 7 = 273.9 n B h 2 . (10) 

The abundance of 4 He is very sensitive to the early expansion rate. Since a non-standard expansion rate (5^1) would result in 
fewer or more neutrons at BBN and, since most neutrons are incorporated into 4 He, the predicted 4 He abundance differs from 
that in SBBN (5=1; AN^ = 0). In contrast, the abundance of 4 He is not very sensitive to the baryon density since, to first order, 
all the neutrons available at BBN are rapidly converted to 4 He. For rjip » 6, N v ~ 3 and, for a primor dial 4 He mass fraction in 
the range 0.23 < Yp < 0.27, to a very good approximation (Kne ller & Stei gman 2004; Steigman 2007), 

Y P = 0.2485 ±0.0006 + 0.0016[(?7io-6)+ 100(5- 1)]. (11) 

In eq. [TT1 the effect of incomplete neutrino decoupling on the 4 He mass fraction is accounted for according to Mangan o et al.l 
(2005) and 5 is related to AN„ (AN^ = N^ - 3.0) by eq. 5. As a result, for a fixed 4 He abundance, a variation in 7710 of ~ ± 0.2 
(corresponding to a ~ 3% uncertainty in the baryon density) is equivalent to an uncertainty in AN, y of ~ ± 0.02. 

In contrast to 4 He, since the primordial abundances of D, 3 He and 7 Li are set by the competition between two body produc- 
tion and destruction rates, they are more sensitive to the baryon density than to the expansion rate. For example, for 7710 ~ 6 
and for a primordial ratio of D to H by n umber, yu = 10 5 (D/H)p, in the range, 2 < yu < 4, to a very good approximation 
dKneller & Steigmarj2004tlSteigman 20071). 



y D = 2.64(1 ±0.03) 



(12) 



L>7io-6(S-l)J 

The effect of incomplete neutrino decoupling on this prediction is at the ~ 0.3% level (Mangan o et al.|[20Q5b . about ten times 
smaller than the overall error estimate above. 



3.1. Observationally-Inferred Primordial Abundances 

Given the monotonic post-BBN evolution of deuterium (as gas is cycled through stars, deuterium is destroyed) and the sig- 
nificant dependence of its predicted BBN abundance on the baryon density (yu? °c ?7io' 6 )> deuterium is the baryometer of choice 
among the light nuclides produced during primordial nucleosynthesis. While observations of D/H in the solar system and the 
local interstellar medium provide a lower limit to the relic deuterium abundance, it is the D/H ratio (by number) measured from 
observations of high redshift, low metallicity QSO absorption line systems whi ch provide an estimat e of i ts primordial abundan ce. 
The weighted m ean of the six, high redshift, low metallicity D/H ratios from lKirkman et all d2003l) and lO'Meara et al] d2006l) is 
(ISteigmanll2007l) 

y DP = 2.68tg. (13) 

Since the post-BBN evolution of 3 He is more complex and model dependent than that of deuterium and, since 3 He is only 
observed in chemically-evolved H 11 regions in the Galaxy and, since the He primordial abundance is only weakly dependent on 
the baryon density (10 5 ( 3 He/H) = y^ oc rf$ 6 ), its role as a baryometer is limited. 
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FIG. 1. — The probability distribution of the baryon density parameter, r)io. The dashed line shows the probability distribution inferred from SBBN (N„ = 3) 
and the adopted primordial abundance of deuterium (see §3). The solid line is the probability distribution of r?io inferred for N„ = 3 from the combination of the 
WMAP-5yr data, small scale CMB data, matter power spectrum data from 2dFGRS and SDSS LRG, SNIa, and the HST Key Project (see §4). 

As for deuterium, the post-BBN evolution of He is monotonic, with Yp increasing along with increasing metallicity. At 
low metallicity, the 4 He abundance should approach its primordial value. As a result, it is the observations of helium and 
hydrogen recombination lines from low-metallicity, extragalactic Hn regions which are most useful in determining Yp. At 
present, corrections for systematic uncertainties (and their unce rtainties) domin ate estimates of the observationally-inferred 4 He 
primordial mass fraction and, especially, of its error. Following Steigman (2007), we adopt for our estimate here, 

Y P = 0.240 ±0.006. (14) 

While the central value of Yp adopted here is low, the conservatively-estimated uncertainty is relatively large (some ten times 
larger than the uncertainty in the BBN-predicted abundance for a fixed baryon density). In this context, it should be noted that 
although very careful studies of the systematic errors in very limited samples of H n regions provide poor estimators of Y P as a 
result of their uncertain and/or mode l-dependent extrap o lation t o zero metallicity, they are of v alue i n providing a r o bust u pper 
bound to Y P . Using t he results of lOlive & Skillman I d2004l) . iFukugida & Kawasaki I (120061) . and iPeimbert et all d2007l) . we 
follow Steigman (2007) in adopting, 

Y P < 0.251 ±0.002. (15) 

Although the BBN-predicted 7 Li relic abundance provides a potentially sensitive baryometer ((Li/H) oc rf , for 7710 > 3), its 
post-BBN evolution is complicated and model-dependent. For these reasons, it is the observations of lithium in the oldest, most 
metal-poor stars in galactic globular clusters and in the halo of the Galaxy which have the potential to provide the best estimate 
of the primordial abundance of 7 Li. The complication associated with this approach is that these oldest galactic stars have had the 
most time to dilute or deplete their lithium surface abundances, leading to the possibility that the observed abundances require 
large, uncertain, and model-dependent corrections in order to infer the prim ordial ab undance of 7 Li. In the absence of corrections 
for depletion, dilution, or gravitational settling, the data of lRvan et al.l d2000l) and lAsplund e t al. (2006) suggest 

[Li] P = 12 + log(Li/H) = 2.1±0.1. (16) 

In contrast, in an attempt to correct for evolution of the surface lithium abundances. iKorn et al.l d2006) use their observations of a 
small, selected sample of stars in the globular cluster NGC6397, along with stellar evolution models which include the effect of 
gravitational settling to infer 

[Li] P = 2.54 ±0.1. (17) 

In the following analysis, the inferred primordial abundances of D and 4 He adopted here are used to estimate 7710 and AN, y . 
Given the inferred best values and uncertainties in these two parameters, the corresponding BBN-predicted abundance of 7 Li can 
be derived and compared to its observationally inferred abundance. 

3.2. BBN Constraints On N u And 7/b 

Since the primordial abundance of deuterium is most sensitive to 7710, while tha t of 4 He is most sensitive to N„, isoabundance 
contours of D /H and Yp in the { 7710, AN„ } plane are very nearly ortho gonal; see lKneller & Steigmani (120041) . The analytic fits 
to BBN from Kneller & Steigman (2004), updated by Steigman (2007), are used in concert with the primordial abundances of 
D and 4 He adopted here to infer the best values, and to constrain the ranges of 7710 and AN^. While these fits do have a limited 
range of applicability, they are, in fact, accurate within their quoted uncertainties for the range of parameter values and observed 
abundances considered here. 

In SBBN, with three species of neutrinos (ANj, = 0), the primordial abundances are only functions of the baryon density, r/s- 
For SBBN, the primordial deuterium abundance adopted in §3.1, yup = 2.68^Jj 25, implies, 

77 I0 (SBBN) = 6.0±0.4 (18) 



FIG. 2. — The 68% and 95% contours in the N„ - 7710 plane derived from a comparison of the observationally-inferred and the BBN-predicted primordial 
abundances of D and 4 He. The shaded region is excluded by the 95% upper bound to the helium abundance in eq. 15 (see, eq. 20). 

This result is in excellent agreement with the independent estimate of 7710 = from the CMB and LSS (discussed below in 

§4). The probability distributions of 7710 inferred from SBBN and from the CMB and LSS are shown in FigureQ] 

For non-standard BBN, with AN, y =/0 (S 1), there is a second free parameter, N„ (or, S). In this case, in addition to yop, the 
4 He abundance Yp is used to constrain the [rjio, } pair. Adopting the D and 4 He abundances from §3.1 (Yp = 0.240 ± 0.006), 
along with the analytic fits in eq. (11) and eq. (12), leads to, 

7710 = 5.7 ±0.4, N„ = 2.4±0.4. (19) 

In Figure |2]are shown the 68% and 95% contours in the N,y-rj W plane which follow from a comparison of the BBN predictions 
with the observationally-inferred primordial abundances of D and 4 He. Notice that while the best fit value of N„ is less than the 
standard-model value of N„ = 3, the standard-model value is consistent with the relic abundances at the ~ 68% level. 

These results are sensitive to the choices of the relic abundances of D and 4 He. We note that if deuterium is ignored and the 
robust upper bound to the 4 He mass fraction, Yp < 0.255 at 95% confidence (eq. 15), is adopted, then eq. 11 provides an upper 
limit to S (Nj,) as a function of the baryon density, 

S< l.lO-O.Olrjio. (20) 
For rjio in the range 5 < 7710 < 7 (see §4), this corresponds to a robust upper bound to N„ ranging from 3.6 to 3.4. 

4. N„ AND r? B FROM THE CMB AND LSS 

The pattern of temperature fluctuations in the cosmic microwave background contain information about the baryon density and 
the radiation density and thus serve as complementary probes of r/s and N u some ~ 10 5 years after BBN. 

The baryon density, parametrized by 773 or Q#h 2 affects the relative amplitudes of the peaks in the CMB temperature power 
spectrum. The ratios of the amplitudes of the odd peaks to the even peaks provide a determination of 773 that is largely uncorrelated 
with N„. 

The radiation density, parametrized by an effective number of neutrino species N„, affects the CMB power spectrum primarily 
through its effect on the epoch of matter-radiation equality. There are substantial differences in amplitudes between those scales 
that enter the horizon during the radiation dominated era and those that enter the horizon later, in the matter dominated era. 
Increasing the radiation content delays matter-radiation equality, bringing it closer to the epoch of recombination, suppressing 
the growth of perturbations. As a result, the redshift of the epoch of matter-radiation equality, z eq , is a fundamental observable 
that can be extracted from the CMB Power Spectrum. z eq is related to the matter and radiation densities by, 

l+Zeq = Pm/PR- (21) 

Since pr depends on N„, z eq is a function of both N„ and £1m^ 2 , leading to a degeneracy between these two parameters. For a flat 
universe, preserving the fit to the CMB power spectrum when N„ increases, requires that f^M and/or Hp increase. As a result of 
this degeneracy, the CMB power spectrum alone imposes only a very weak constraint on N„ dCrottv et al JI2003I : lPierpaolj||2003b 
Barger et al. 20031: lHanne stad 2003; Ichikawa et al. 2006). Inclusion of additional, independent constraints on these parameters 
are needed to break the degeneracy between N„ and Q M h 2 - 

Besides affecting the epoch of matter-radiation equality, relativistic neutrinos leave a distinctive signature on the CMB power 
spectrum due to their free streaming at speeds exceeding the sound speed of the photon-baryon fluid. This free streaming creates 
neutrino anisotropic stresses generating a phase shift of the CMB acoustic oscillations in both temperature and polarization. This 
phase shift is unique and, for adiabatic initial conditions, cannot be generated by non-relativistic matter. In principle, this e ffect 
can be used to break the degeneracy between N„ and tlyih 2 , leading to tighter constraints on N„ dBashinskv & Seliakll2004l) . 
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FIG. 3. — The top panel shows the CMB power spectrum for the best fit models with N„ fixed at N„ = 1 (solid black), N„ = 3 (dashed red), and N„ = 5 
(dot-dashed blue) illustrating its insensitivity to N„ in the absence of an independent constraint on Q.yih 1 . The bottom panel shows the matter power spectra for 
the same set of parameter values, illustrating its sensitivity to N u . 

Alternatively, since the luminosity distances of type la supernovae (SNIa) provide a constraint on a combination of 17m and fl\ 
complementary to that from the assumption of flatness, they are of value in restricting the allowed values of SIm- In concert with 
a bound on Hq, this, too, helps to break the degeneracy between N„ and Vly[h 2 . 

Another way to break the degeneracy between N„ and VLy[h 2 is to use measurements of the matter power spectrum in com- 
bination with the CMB power spectrum. To preserve a fit to the CMB power spectrum, an increase in N„ requires that Sl^h 2 
increase in order that the redshift of matter-radiation equality remain unchanged. The turnover scale in the matter power spec- 
trum is set by this connection between N v and SIm^ 2 - Since the baryon density is constrained by the CMB power spectrum, 
independently of N„, increasing the radiation density (N„ > 3) requires a higher dark matter density in order to preserve z eq (in 
a flat universe, f^M + ^A = 1)- Between the epoch of matter-radiation equality and recombination the density contrast in the cold 
dark matter grows unimpeded, while the baryon density contrast cannot grow. Consequently, increasing N„ and Vlyji 2 , increases 
the amplitude of the matter power spectrum on scales smaller than the turnover scale corresponding to the size of the horizon at 
z eq . Data from galaxy redshift surveys can be used to infer the matter power spectrum, thereby constraining fl M h 2 and N„ . This 
effect may be seen in Figure [3] which shows that for nearly indistinguishable CMB power spectra, different values of N u yield 
distinguishable matter power spectra. The upper panel of Figure[3]shows that by making suitable adjustments to the other cosmo- 
logical parameters, specifically the matter density and the spectral index, CMB power spectra which are nearly degenerate up to 
the third peak of the power spectrum can be produced using very different values of N„. However, as the bottom panel of Figure 
[3]illustrates, these models produce matter power spectra with different shapes, demonstrating that the matter power spectrum can 
be used to help break the degeneracy between N„ and UmIi 2 . 

4. 1 . Analysis And Datasets 
Before presenting our results, we describe the analysis and the datasets employed. 

For our analysis we assume a flat, CDM cosmology with a cosmological constant, A, and three flavors of active neutrinos with 
negligible masses. Our cosmological model is parametrized by seven parameters: 

p = {n B h 2 ,n M h 2 ,h,T,n s ,As,K v }. (22) 

The contents of the universe are described by the baryon density, 51b/j 2 and the matter (baryonic plus cold dark matter) density, 
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Cluh ■ Since a flat cosmology is assumed, the dark energy density and the matter density are related by fl\ = 1 - £1m- The expan- 
sion rate of the universe is described by the reduced Hubble parameter, h (Hq = lOOh kms _I Mpc _1 ). Instantaneous reionization 
of the universe is assumed with optical depth to last scattering, r. As is the amplitude of scalar perturbations and n$ is the scalar 

spectral index. 

The Code for Anisotropics in the Microwave Background (CAMB) dLewis et al.ll2000l) is used to compute the CMB power 
spectrum for a fixed set of cosmological parameters. For a given dataset, our degree of belief in a set of cosmological parameters 
{p} is quantified by the posterior probability distribution, 

P(p | data) oc L(data \ p) Tl(p) . (23) 

The likelihood L(data|p) quantifies the agreement between the data and the set of parameter values {p}. Il(p) represents the prior 
on cosmological parameter values before the data are considered. 

Markov Chain Monte Carlo meth ods are used to explore t he multi-dimensional likelihood surface. We use the publicly available 
COSMOMC code for our analysis dLewis & Bridle 1 12002). Flat priors are adopted for all parameters, along with a prior on the 
age of the universe, to, of to > 10 Gyr. 

The mode of the marginalized posterior probability distribution is used as a point estimate and the minimum credible interval 
as an estimate of the uncertainty. The minimum credible interval selects the region of the parameter space around the mode, 9, 
that contains the appropriate fraction of the volume (e.g., 68%, 95%) of the posterior probability distribution, while minimizing 
9-6. The minimum credible interval selects the region of the parameter space with the highest probability densities 3 . 

Our primary d ataset is the CMB data from the Wilkin son Microwave Anisotropy Probe (WMAP) accumulated from five years 
of observations (Hinshaw et al. 2008; Nolta et al. 2008). For the WMAP data, likelihoods are computed using the code provided 
on the LAMBDA webpage 4 . There are a number of ground and balloon based CMB experiments whose high angular resolution 
probe smaller scales than those probed by WMAP. These are more sensitive to the higher order acoustic oscillations beyond the 
third peak i n the CMB power spe ctrum. In particular, the 2008 results from the Arcminute Cosmology Bolometer Array Receiver 
(ACBAR) (Reic hardt et al.l2008l) impose the strongest constraints at present on the CMB p ower spectrum at small angular scales. 
In ou r analysis, data from the follow ing C MB experiments are used: B OOMERANG dPiacentini et alj 120061: iMontrov et all 
|2006|), ACBAR ^^S51]^B i-^L^^^ ^51]|^^ IMaim et al.ll2003l) . VSA dDickinson et alJl2004 . MAXIMA 
(iHanany et al.ll2000h and DASI (lHalverson et al.ll2001l) . 

To help break the degeneracy between N„ and ^mIt 2 , w e adopt a Gaussian p rior on the Hubble parameter, Ho = 72 ± 
8 kms _I Mpc _1 from the Hubble Space Telescope Key Project dFreedman et al.ll2~001l) . 

Since the luminosity distances of type la supernovae (SNIa) provide a constraint on a combination of f?M and 51a, they are of 
value in restricting the allowed values of Om, helping to break the degeneracy between N„ and Vluh 2 . W e use the luminosity 
distances for 115 type la supernovae measured by the Supernova Legacy Survey (SNLS) ( Ast ier et al.l2006l) . For each supernova, 
the observed luminosity distance is compared to that predicted for a given set of cosmological parameters. 

We use the matter powe r spectrum inferred from galaxy re dshift surveys such as the Sloan Digital Sky Survey (SDSS) 
(Teg marketafl 2004, 2006) and the 2dFGRS (ICole etalJ 120051) . In using the matter power spectrum to infer cosmological 
parameters, it has become clear that the constraint on fijyi is sensitive to the choice of length scales on which the power spec- 
trum is measured. The power spectrum on smaller scales prefers higher v alues of f^M, provided that those scale s are correctly 
described by linear perturbation growth and scale independ ent galaxy bias dCole et alj|2005t iPercival et al.ll2007|). For example . 
using the SPSS LRG power spectrum dTegmark et al.ll2006h in combination with the 3 year WMAP dataset dSpergel et alj 2007). 

(2Q08) find a disagreement between the matter density inferred on scales with k < O.lh Mpc" 1 and k < 0.2h Mpc" 1 . 
(2007) have shown that the constraint on N„ obtained from the matter power spectrum is affected similarly, result- 



Dunkl ev et al 



Hamann et al 



ing from the correlation between S1m and N, y . Therefore, here we only use matter power spectrum data on scales that are likely to 
be safely linear. We truncate the matter power spectrum at k = 0.07 h Mpc" 1 , keeping data only on scales with k < 0.07/z Mpc" 1 . 
We assume that galaxy bias is constant and scale independent for these scales. 

4.2. N v And The Lyman-a Forest 

Measurements of the flux power spec trum of the Lyman- a forest in QSO ab sorption spectra can be used to reconstruct the 
matter power spectrum on small scales dCroft et alJ fl998: McD onald et al.ll2000|). Observations of the SDSS Ly-a flux power 
spectrum have been used to constrain the linear matter power spectrum at z ~ 3 dMcDonald et"ai1l2005l) . When these constraints 
are comb ined with CMB po wer spectrum data, they favor considerably higher val ues of N„ than those we find here (§5). For 
example. [Seljak et all (120061) found N„ = 5.2 and a 95% range 3.4 < N„ < 7.3 and lHamann et al l d2007l) found N„ = 6.4 and a 
95% range 3.2 < Nj, < 11, the difference being accounted for by their different combinations of data sets. Each of these excludes 
the standard model value of N„ = 3 at more than 95% confidence. 

As discussed earlier and shown in the lower panel of Figure [3] the principal effect of an increase in N v is to increase the 
amplitude of the matter power spectrum on scales smaller than those corresponding to the horizon at matter-radiation equality, 
Zeq- The ACDM fits to the Lyman-a forest data prefer higher amplitudes of density fl uctuations on small scales compared to those 
expected from measurements of the WMAP power spectrum dViel & Haehneltll2006l) . favoring higher values of N„. This effect is 
seen in Figure [3] which shows that for nearly indistinguishable CMB power spectra, different values of N„ yield distinguishable 
matter power spectra. To preserve z eq and the fit to the CMB, the model with N„ = 1 has a lower value of Q,uh 2 and, therefore, 

3 The commonly used GETDIST analysis packa ge uses the mean o f the marginalized posterior probability distribution as a point estimate and gives uncertainty 
estimates based on the central credible interval (Hamann et al. 2007). These estimates are identical for Gaussian probability distributions but differ significantly 
for non-Gaussian distributions, particularly for asymmetric probability distributions. 

4 http://lambda.gsfc.nasa.gov 
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FIG. 4. — The 68% and 95% contours in the N„ - Hm/i 2 plane inferred from the combination of the WMAP-5yr data, small scale CMB data, luminosity 
distances of SNIa and the HST Key Project prior on Hq. The dashed line shows the locus of points corresponding to the same value of Zeq(= 3144), illustrating 
the degeneracy between these two parameters. As the contours reveal, this degeneracy may be broken if complementary data is used to constrain SIm/i 2 . 

the matter power spectrum for that model has a lower amplitude on scales smaller than the scale corresponding to the horizon 
at the epoch of matter-radiation equality. Conversely, the model with N„ = 5 has a higher value of Om/i and the matter power 
spectrum for that model has a higher amplitude on scales smaller than the scale corresponding to the horizon at z eq . 

Before reaching any conclusions about N„ based on the Lyman-a forest data, it is worth noting that assumptions about the 
thermal st ate of the IGM play an important role in reconstructing the matter power spectrum from the Lyman- q forest flux power 
spectrum. Bolt on et"aLl (120071) compare measurements of the Lyman-a forest flux probability distribution by iKim et al.l (120071) 
to hydrodynamic simulations of the Lyman - q fore st, finding evidence for an inverted temperature-density relation for the low 
density intergalactic medium. iBolton et al.l (120071) suggest that He n reionization could be a possible physical mechanism for 
achieving an inverted temperature-density relation. Such an inversion would result in a smaller amplitude of the matter power 
spectrum for a given observed flux power spectrum, thereb y alleviating the tensio n with the other data sets which tended to drive 
N„ to high values. However, in their power spectrum fits lMcDonald et all (120051) marginalize over equation of state parameters, 
so this explanation of the tension may not be entirely satisfactory. Future studies of the Lyman-a forest may weaken or strengthen 
the evidence for a discrepancy. For these reasons, we do not use data from the Lyman-a forest in our analysis. 

4.3. CMB And LSS Constraints On N v And ?7b 

According to iDunklev et"aTI (l2008h the WMAP 5 year data only impose a lower limit on of > 2.3 but, according to 
iKomatsu et al.l (120081) . do not lead to an upper limit due to the degeneracy between N„ and Vly[h 2 . Inclusion of data from small 
scale CMB experiments do not break this degeneracy. Figure [4] illustrates this degeneracy, as well as how it may be broken by 
non-CMB constraints on Vly[h 2 . The dashed line in Fig.|4]is the locus of points with constant z eq = 3144. Figure|4]shows the joint 
probability distributions of N„ and Q,yji 2 inferred from the WMAP 5yr data and small scale CMB experiments, supplemented 
by independent data from measurements of SNIa luminosity distances and the HST Key Project prior on Hq are used to bound 
Vly[h 2 . The range of N„ is now limited. 

Our constraints on N, y and 7/10 from various CMB and LSS datasets and combinations of them are summarized in Table 1 and 
in Figures 5-7. 

Using the WMAP 5 year data in combination with data from other CMB experiments (ACBAR, BOOMERANG, CBI, DASI, 
MAXIMA and VSA), along with the HST Key Project prior on Hq and luminosity distance measurements of type la supernovae, 
we find central values Nj, = 2.9 and r] W = 6.2, along with 68% (95%) ranges of 2.0 < N^ < 4.1 (1.3 < N„ < 5.4) and 6.0 < r] W < 
6.3 (5.9 < 7710 < 6.4). Figure [5] shows th e joint probability distribution of N^ and 7710 for this combination of CMB datasets. 

Adding the 2dFGRS power spectrum (ICole et al .1120051) to the CMB power spectrum from the WMAP 5 year dataset, along 
with ground based CMB experiments (see above), the HST prior on Hq and, luminosity distance measurements from SNIa, we 
find central values N„ = 3.0 and 7710 = 6.1, along with 68% (95%) ranges 2.1 <N„< 4.2 (1.3 <N„< 5.2) and 6.0 < 7710 < 6.3 
(5.9 < 7710 < 6.4) re spectively; see Table 1. Replacing the 2dFGRS power spectrum with the SDSS LRG power spectrum 
dTegmark et al.ll2006l) . we obtain very similar resu lts: N„ = 2.8, 2.1 < N „ < 3.9 (1.5 < N„ < 5.2) and, rj l0 = 6.2, 6.1 < 7710 < 6.3 
(5.9 < 7710 < 6.5). In contrast, the SDSS DR2 (Tegmarket al. 2004) power spectrum favours considerably higher values of 
Ny compared to those inferred from the 2dFGRS and the SDSS (LRG) power spectra. 

Our best constraints on N„ and 7710 from the CMB and LSS data are derived by combining the following datasets - WMAP 
5-year data, ground and balloon based CMB experiments (BOOMERANG, ACBAR, CBI, VSA, MAXIMA and DASI), the 
HST prior on Ha, the SNIa luminosity distance measurements and, the matter power spectrum on large scales inferred from the 
2dFGRS and SDSS (LRG) data. Constraints on the matter power spectrum from the Lyman-a forest are not included for the 
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FIG. 5. — The 68% and 95% contours in the N„ — t][q plane inferred from the combination of the WMAP-5yr data, small scale CMB data, SNIa luminosity 
distances, and the HST Key Project prior on Hq (see the text and Table 1). 




FIG. 6. — The 68% and 95% contours in the N„ —7710 plane derived using the WMAP 5-year data, small scale CMB data, SNIa and the HST Key Project prior 
on Hq along with matter power spectrum data from 2dFGRS and SDSS LRG(see the text and Table 1). 



reasons discussed above in 34.21 Using these datasets, we obtain (see Table 1): 



ma 



:6.1 



+0.2+0.3 
-0.1-0.2- 



M _ ?q +l. 0+2.0 
IN!/ - -i-y_o.g_i.4- 



(24) 
(25) 

At present the combined CMB and LSS data provide the best baryometer, determining the baryon density to better than 3%, 
but only a relatively weak chronometer, still allowing a large range in S (0.87 < S < 1.14 at 95% confidence). Within their 
uncertainties, the CMB/LSS data, which probe the Universe at > 400,000 years, are consistent with BBN, which provides a 
window on the Universe at < 20 minutes. 

5. COMPARING THE BBN AND CMB/LSS CONSTRAINTS 

Using the ranges of 7710 and N„ allowed by the CMB and LSS data, the BBN-predicted primordial abundances of 4 He, D and 7 Li 
may be inferred. Figure|8]compares these constraints to the observationally-inferred primordial abundances adopted in §3.1. The 
CMB/LSS -inferred BBN abundances of D and 4 He are in excellent agreement, within the errors, with the observationally-inferred 
relic abundances. For the central values of 7710 = 6.14 and N„ = 2.9, the BBN-predicted deuterium abundance of yop = 2.54 is, 
within the errors, in agreement with its observationally-inferred primordial value of yop = 2.68. For 4 He, the BBN-predicted 
mass fraction is Yp = 0.247, slightly high compared to the central value of the primordial abundance adopted in §3.1, Yp = 
0.240, but within 1.2cr of it and, completely consistent with the evolution-model independent upper bound presented in eq. 15, 
Y P < 0.251 ±0.002. 

For 7 Li the BBN-predicted best fit from the CMB/LSS d ata is [Lib = 2.66, c onsiderably higher tha n the value ([Li] P = 2. 1 ±0. 1) 
determined from observations of metal-poor halo stars dRvan et alj d2000l) . lAsplund et al.l d2006l)) w ithout any correction for 
depletion, destruction, or gravitational settling. If, however, the correction proposed bv lKorn et al.l ( 2006) is applied, the predicted 
and observed 7 Li abundances may, perhaps, be reconciled, as may be seen from the lower panel of Fig. [8] It remains an open 
question whether this lithium problem is best resolved by a better understanding of stellar physics or, if it is providing a hint of 
new physics beyond the standard model. 
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TABLE 1 

N„ AND 7710 FROM DIFFERENT DATASETS 



Dataset 




N„ 


J?10 


BBN (Y P & yr> P ) 


2 


,1+0.4+0.9 
-0.4-0.8 


c 7+0.4+0.8 
J ' -0.4-0.9 


WMAP(lyr)+HST 


2 


0+4.5+5.5 
—1.8-1.9 


ft 1+0.5+1.8 
u-1 -0.3-0.7 


WMAP(3yr)+HST+SN 


2 


Q+2. 1+4.0 

2-2.1 


6 0fi 4O.23+0.42 
-0. 20— 0.39 


WMAP(3yr)+CMB(07)+HST+SN 


2 


c+1.7+3.8 
-1.2-1.7 


ft , -1+0.21+0.39 
u - 1J -0.18-0.36 


2dFGRS+WMAP(3yr)+CMB(07)+HST+SN 


2 


Q+l. 4+2.3 
-1.2-2.3 


ft 1 1+0.20+0.38 
U-11 -0. 15-0.32 


SDSS(DR2)+WMAP(3yr)+CMB(07)+HST+SN 


3 


7+1.6+3.3 
-1.4-2.4 


ft ic+0. 15+0.35 
u - 1J -0.20-0.39 


SDSS(LRG)+WMAP(3yr)+CMB(07)+HST+SN 


2 


n +l. 2+2.6 
"-1.2-1.8 


ft n+0.19+0.36 
°- lz -0.16-0.34 


WMAP(5yr)+HST+SN 


3 


q+2.0+4.5 
-1.2-2.4 


ft 1 1+0.16+0.33 
u - 11 -0.16-0.33 


WMAP(5yr)+CMB+HST+SN 


2 


n+l. 2+2.5 
'-0.9-1.6 


ft 16 +0.14+0.25 
u - lu -0.16-0.30 


2dFGRS+WMAP(5yr)+CMB+HST+SN 


3 


n+i. 2+2.2 
u -0.9-1.7 


ft 1^+0.16+0.27 
a - 14, -0. 14-0.25 


SDSS(LRG)+WMAP(5yr)+CMB+HST+SN 


2 


c+1.1+2.4 
-0.7-1.3 


6 jft+0.14+0.27 
u ' lu -0.14-0.27 


SDSS(LRG)+2dFGRS+WMAP(5yr)+CMB+HST+SN 


2 


n+l. 0+2.0 
-0.8-1.4 


ft 1 4+0. 16+0.30 
u - 1H, -0.11-0.25 


BBN+SDSS(LRG)+2dFGRS+WMAP(5yr)+CMB+HST+SN 


2 


c+0.4+0.7 
J -0.4-0.7 


ft 1 1+0.12+0.26 
°- 11 -0. 13-0.27 



NOTE. — Best fits and 68% and 95% confidence intervals for N„ and rjio from our prin- 
cipal datasets. WMAP refers to the CMB power spectrum data from the WMAP experi- 
ment and CMB to the data from ACBAR+ BOOM+CBI+VSA+MAXIMA+DASI. CMB(07) 
uses the 2007 A CBAR dataset Kuo et al.ll2007l) while CMB uses the 2008 ACBAR dataset 
iReichardt et al. 2008). HST refers to the prior on H from the HST Key Project. SN stands 
for the SNIa luminosity distance measurements. SDSS (LRG) and 2dF refer to the respective 
LSS matter power spectra, truncated at k = 0.07/1 Mpc -1 . 




FIG. 7.— 

On the left (red, dashed), the probability distribution of ti v inferred from the combination of the WMAP 5-year data, small scale CMB data, SNIa and the HST 
Key Project prior on Ho and matter power spectrum data from 2dFGRS and SDSS LRG. The solid blue curve is the BBN (D plus 4 He) distribution. On the right 
(same line types and colors), are the probability distributions of »jio using the same data sets. 

6. CONCLUSIONS 

While our CMB+LSS constraints on rjn and are consistent with most previous analyses (IjCrotty et al.l 12003b iPierpaolil 
20031; iHannestadl 120031: iBarger et all 1207)1 ISeliak et al.l|2006t llchikawa etaT1l2006t ISpergel et al] 120071; iMangano et al J 120071 
Hamann et alj|2007l ; iDunklev et al. 2008; K omatsu et al.ll2008r) . they are tighter because we have used more and/or more recent 
data. However, as discussed above in 34.21 analyses that included the Lyman-a forest data generally find higher values of N^. 

Until the advent of WMAP and the other ground- and balloon-based CMB experiments, BBN provided the best baryometer 
(mainly from deuterium) and chronometer (mainly from helium-4). As may be seen from Table 1, while the WMAP first year 
data provided a competitive baryometer, it offered a relatively poor chronometer. This improves with the WMAP 3 and 5 year 
data, especially when they are combined with the other CMB and LSS data. These now lead to a determination of the baryon 
density at the ~ 2-3% level, a factor of > 2 better than that from BBN. However, although the CMB/LSS constraint on N^ has 
improved significantly and, it is consistent with that from BBN, it remains weaker than the BBN constraint by a factor of ~ 2. 
For some time now BBN has clearly established at high confidence that N, y > 1 w hen the Universe w as ~ 20 minutes old. For 
example, using a sli ghtly different est imate of the primordial helium mass fraction. [Barger et al.l (120031) found N u > 1.7 at 95% 
confidence (see also Steigman (2006a)). The more recent WMAP 5 year data, combined with the recent ACBAR and other CMB 
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Fig. 8.— 

The solid black curves show the probability distributions for the primordial abundances of D, 4 He and 7 Li derived from the values of rj^ and N„ inferred from the 
CMB, LSS, SNIa and the HST prior and matter power spectrum data from 2dFGRS and SDSS LRG. The blue dashed curves show the probability distributions 
for the observatio nally-inferred primordial abundances of D, 4 He and 7 Li; see §3.1. The red, dot-dashed curve in the far right panel is the 7 Li abundance from 
IKometaljEOOa . 
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Fig. 9. 



(Left) In blue (solid), the 68% and 95% contours in the N„ - 7710 plane derived from a comparison of the observationally-inferred and BBN-predicted primordial 
abundances of D and 4 He (see Figure 2). In red (dashed), the 68% and 95% contours derived from the combined WMAP 5-year data, small scale CMB data, 
SNIa, and the HST Key Project prior on Hq along with matter power spectrum data from 2dFGRS and SDSS LRG (see the text and Table 1). (Right) The 68% 
and 95% joint BBN-CMB-LSS contours in the N„ -7)10 plane. 

and LSS datasets, the CMB/LSS now confirm that N„ > 1 (or, N„ > 2 dDunklev et alj |2008)) when the Universe was > 400 kyr 
old. 

As may be seen from Table 1 and from Figures [7] [8] and from the left-hand panel of Figure[9] BBN and the CMB, which probe 
physics at widely separated epochs in the evolution of the Universe, are in excellent agreement. This permits constraints on any 
differences in physics between BBN and recombination and/or the present epoch. For example, since baryons are conserved, t/b 
relates the number of thermalized black body photons in a comoving volume at different epochs (see figure |7), constraining any 
post-BBN entropy production. Our results from the CMB/LSS and from BBN imply, 



jyCMB 

-^- = 0.92 ±0.07. 

N BBN 



(26) 



This ratio is consistent with 1 at ^ler and places an interesting upper bound on any post-BBN entropy production. 

Alternatively, late decay ing particles could pro duce relativistic particles (radiation), but not necessarily thermalized black body 
photons (see, for example, Ichikaw a et alj d2007l) ). Deviations from the standard model radiation density can be parametrized by 
the ratio of the radiation density, p' K to the standard model radiation density, pr. In the post-<? ± annihilation universe (see eq. 9), 



R = S 2 = ^ = 1+0.134 AN„ 

PR 



(27) 



Comparing this ratio, at BBN and at recombination (see figure |7), any post-BBN production of relativistic particles can be 
constrained. 



^CMB 
^BBN 



- 1 07 +016 



(28) 



This ratio, too, is consistent with 1 within ler, placing an upper bound on post-BBN production of relativistic particles. 
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Although the non-standard expansion rate has been parametrized in terms of an equivalent number of additional species of 
neutrinos, we have emphasized that a non-standard expansion rate need not be related to extra (or fewer) neutrinos. For example, 
deviations from the standard expansion rate cou ld occur if the value of the early-Universe gravitational constant, Gn were different 
in the from its present, locally-measured value (Yang et al.lll979l : lBoesgaard & Steig man 1985t lAccetta et alll990tlCvburt et al.l 
2005). For the standard radiation density with three species of light, active neutrinos, the constraint on the expansion rate can be 
used to constrain variations in the gravitational constant, Gn- From BBN, 

qBBN 

S 2 = ^— = 0.91 ±0.07 (29) 
Gn 

and at the epoch of the recombination, 

5 2 = ^ = 099! o.n (30) 

consistent with no variation in G at the ~lcr level. 

The agreement between rjs and N„ evaluated from BBN (~ 20 minutes) and from the CMB/LSS (> 400 kyr) is shown in the 
left hand panel of Figure|9] As Figure|9]illustrates, BBN and the CMB/LSS, which probe the Universe at widely separated epochs 
in its evolution, are completely consistent. As already noted, at present the CMB is a better baryometer while BBN remains a 
better chronometer. Since these independent constraints from the CMB/LSS and BBN are in very good agreement, we may 
combine them to obtain the joint fit in Table 1 and shown in the right hand panel of Figure|9] We note that while the best-fit value 
of N^ is less than 3, this is not statistically significant since the results are consistent with the standard model of three species of 
active neutrinos at the ~lcr level. 

Of course, our BBN results are sensitive to the relic abundances we have adopte d. For comparison, w e have repeated our 
analysis for an alternate set of primordial abundances. For deuterium, we adopted the O'Me ara et al.l d2006) results based on the 
weighted mean of the log(yo) values, yo = 2.84^ |j and, for 4 He we chose, somewhat arbitrarily, Yp = 0.247 ±0.004. For this 
alternate abundance set the BBN-predicted value of t]b is virtually unchanged from our previous result, 7710 = 5.7 ± 0.3, while 
Nj, = 2.9 ±0.3 is much closer to the standard model expectation. As a result, while the constraint on entropy production (eq. 26) 
is unchanged, Rcmb/Rbbn = l-00^ {f and G^ BN /G N = 0.99^;^. When the alternate BBN constraints are combined with those 
from the CMB and LSS, we find r) W = 6.09^g;}| and N„ = 2.9 ± 0.3. 

Very recent observations of deuterium in the high redshift, low metallicity Damped Lyman-a absorber bv lPettini et al.l (2008) 
lead to a deuterium abundance very close to the mean of the previous six abundances used in this paper. As a result, the change 
in yo is very small (yo = 2.70 rather than yo = 2.68 adopted in this paper). Consequently, the change in the parameters inferred 
from it (N„ and 7710) is also very small. 

Future CM B experiments will improv e the constraint on N v by measuring the neutrino anisotropic stress more accurately. 
According to iBashinsky & Seljakl (120041) . PLANCK should determine N, y to an accuracy of <j{N v ) ~ 0.24 and CMBPOL, a 
satellite based polarization experiment, might improve it further to (t(N u ) ~ 0.09, independent of the BBN constraints. In this 
context, we note that such tight constraints on N v will be sensitive to the value of the 4 He abundance adopted in the CMB analysis. 
To achieve these projected accuracies, it will no longer be sufficient to fix Yp in advance. Rather, Yp should be solved for in 
concert with the other cosmological parameters. To this end, we point out that for 7710 m 6, N„ m 3 (su ggested by our results), a 
very good, simple approximation to Yp is provided by eq. 1 1 dKneller & S teigman 2004; Steigman 20Q7|). 
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